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I.  ^INTRODUCTION 

The  influence  of  reagent  vibrational  energy  on  the  rate  of  chemi- 


1 -8 


cal  reactions  is  the  subject  of  a great  deal  of  interest. _ It  is  of 


3 ~ 6 

special  importance  to  laser  enhanced  reactions, atmospheric  reac- 
tions under  disturbed  conditions  and  combustion  processes  in  general. 

4-8 

Recent  studiesv__^have  demonstrated  rate  accelerations  of  several 
orders  of  magnitude  when  one  of  the  collision  partners  is  vibrationally 
excited.  The  reaction  O + H^^  OH  + H has  recently  received 

attention  because  of  its  importance  in  modeling  IR  radiation  in  certain 
rocket  plume  s.~\Birely  et  al.1  have  set  an  experimental  upper  bound  on 
the  rate  for  this  /reaction,  and  Light9  is  currently  measuring  the  actual 

rat6'  l Are  report*^ 

We  report  l^ere  the  results  of  an  extensive  classical  trajectory 
study  of  the  reaction  of  atomic  oxygen  with  both  H^fv^)  and  h£(v=1) 
using  an  LEPS  potential  energy  surface.  In  addition  to  determining 
the  rates  for  these  reactions,  we  have  calculated  the  branching  ratio  i£  calcula” 
for  the  production  of  OH(v'  = l)  and  OH(v'=0)  when  ri^(v=l)  is  the  reactant. 
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II.  POTENTIAL  ENERGY  SURFACE 

The  potential  energy  surface  used  in  these  calculations  was  an 
LEPS  function10  with  a single  adjustable  Sato  parameter.  The  Sato 
parameter  was  adjusted  by  trial  and  error  so  that  the  computed  rate 
constant  of  the  reaction  O + Hg — ►OH  + H would  approximately  equal 
the  experimentally  observed  rate  at  a particular  temperature.  We  are 

O 

primarily  interested  in  low  temperature  results.  At  320  K,  Campbell 

11  - 1 7 

and  Thrush  measured  the  rate  constant  to  be  k = (2.0  ± 0.  16)  X 10 


cm^  molecule  * sec 


In  a very  recent  set  of  low  temperature 

12 


measurements  from  347  to  742  K,  Dubinsky  and  McKenney  state  that 

their  results  extrapolate  to  within  11  percent  of  the  Campbell  and  Thrush 

value.  The  actual  extrapolated  value,  using  their  Arrhenius  parameters, 
-17  3 -1  -1 

is  k = 1.77  X 10  cm  molecule  sec  . The  value  of  k calculated 

13  -17  3 

from  the  Leeds  formula  at  this  temperature  is  k = 0.8  x 10  cm 
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molecule-1  sec-1.  This  is  probably  too  low,  while  the  Campbell  and 
Thrush  point  may  be  somewhat  high.  The  value  of  the  Sato  parameter, 
adjusted  to  make  the  calculated  k fall  somewhere  between  these  limits, 
is  A = 0.0885.  The  calculated  rate  constant,  using  this  value,  is  k = 1.4 
X 10  cmJ  molecule  sec  . This  rate  was  computed  with  in  the 
v = 0 state.  This  is  valid,  since  contributions  of  the  v = 1 state  to  the 

O 

rate  is  negligible  for  temperatures  less  than  1000  K. 

A complete  list  of  the  LEPS  potential  parameters  is  given  in 

Table  I.  A contour  plot  of  this  potential  for  the  collinear  configuration 

is  shown  in  Figure  1.  The  minimum  energy  path  profile  is  displayed 

in  Figure  2.  The  reaction  barrier  along  this  path  is  12.5  kcal/mole 

with  respect  to  the  reactants  channel.  Also  indicated  in  Figure  2 are 

th. 

the  energies  for  the  zero — and  first  vibrational  levels  of  both  products 
and  reactants.  The  barrier  height  is  significantly  above  the  experi- 
mental activation  of  8.  9 kcal/mole.  Of  course  by  adjusting  the  surface 
to  reproduce  the  experimental  rate,  we  have,  to  some  extent,  incorpor- 
ated the  effects  of  zero  point  energy  differences  and  tunneling  into  our 
semiempirical  surface. 

III.  COMPUTATIONAL  PROCEDURE 

A slightly  modified  version  of  Muckerman's  classical  trajectory 
program  CLASTR**  was  used  to  compute  the  state -to -state  reactive 
cross  sections.  The  modification  we  made  consisted  of  changing  the 
method  of  choosing  the  impact  parameters.  Following  a suggestion  of 
Porter,17  the  values  of  the  impact  parameter,  b,  were  restricted  to 
the  discrete  set  that  corresponds  to  integral  values  of  the  orbital  angular 
momentum  quantum  number  I.  Accordingly, 


b,  = 1i(2nE)-1/2  (I  + 1/2), 


where  p is  the  reduced  mass  of  the  system  and  E is  the  initial  relative 
translational  energy.  The  state -to -state  reactive  cross  section  is  then 
given  by 


Table  I 


Potential  Parameters  Used  in  the  LEPS  Surface  for  the  0-Ha~H^  System 


Parameter 

G-Ha 

H -H, 
a b 

0-Hb 

i 

0—1 
e (a  ) 

2.294 

1.  942 

2.  294 

D (kcal/mole) 

106.  6 

109.4 

106.  6 

ro  <A> 

0. 9706 

0. 7417 

0.  9706 

A 

0.0885 

0.0885 

0.  0885 

HH  DISTANCE  IAI 


REACTION  COORDINATE  IAI 


Figure  2.  Schematic  diagram  of  the  energy  along  the  minimum  energy 
path  on  the  potential  surface  shown  in  Figure  1.  Also  indi- 
cated are  the  energy  levels  of  the  zero-  and  first  vibrational 
levels  of  both  products  and  reactants. 
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‘max 


7q(v,  V;  E)  = (nli  /2p,E)  £ (21  + 1)  P(v,  J,  v1,  b^;  E) 

l = 0 


(2) 


where  v and  J are  the  initial  vibration  and  rotation  quantum  numbers 

and  v'  is  the  final  vibration  quantum  number.  (The  final  rotation 

states  were  not  distinguished  in  these  calculations.  ) The  probability, 

P(v,  J,  v' , bp  E),  is  the  ratio  of  the  number  of  reactive  trajectories 

that  terminate  in  the  vibration  state  v1  to  the  number  that  started  in 

the  initial  state  specified  by  v,  J,  and  bp  The  f values  are  varied  in 

unit  steps  from  0 to  f , while  the  remaining  initial  conditions  were 

randomly  selected  by  a Monte  Carlo  procedure.  The  most  convenient 

method  for  determining  fmax  is  to  select  a maximum  impact  parameter 

b and  then  compute  t as  the  nearest  integral  solution  of  Eq.  (1). 
max  max 

The  value  of  b is  determined  by  trajectory  calculations  to  be  the 
smallest  impact  parameter  for  which  no  reaction  occurs. 

Assuming  the  translation  and  rotation  degrees  of  freedom  are  in 
thermal  equilibrium  at  a temperature  T,  the  rate  constant  from  a 
specific  vibrational  state  of  the  reactant  is  computed  from  the  equation 


k(v,  T)  = Q'1  f(T)  |£gj(2J  + 1)  exp(Ev  J/kT) 


00 


X (8kT/TT^)1/2  (l/kT)y cr(v,  J;  E)  exp(-E/kT)  EdE 

where  E T is  the  vibration- rotation  energy  of  the  reactant  molecule, 

gj  is  the  statistical  weight  of  the  rotation  state  J (in  the  case  of  H^, 

g = 1 for  even  J and  g = 3 for  odd  J),  Q is  the  rotational  partition 
J J 

function,  f(T)  is  the  probability  that  the  system  is  initially  on  an 

electronic  surface  on  which  the  reaction  can  occur  (multiple  surface 
18  19 

coefficient  ' ) and  a(v,  J;  E)  is  the  reaction  cross  section  summed 

over  final  vibration  states 


(3) 
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o(v,  J;  E)  = £ an(v’  J»  v';  E)* 
v' 


(4) 


In  order  to  carry  out  the  Integration  in  Eq.  (3)  and  also  to  help  average 
out  statistical  errors,  the  computed  cross  section  points  are  fitted  by 
least  squares  to  the  function 


ct(v,  J;  E) 


0 E £ E^fv,  J) 

4 

£ A (v,  J)  [E'n  - E“"(v,  J)]  E > E^v,  J). 
n=0  J 


(5) 


This  is  substituted  into  Eq.  (3)  which  is  then  evaluated  numerically. 

The  threshold  energy,  E^fv,  J)»  is  determined  by  fitting  a straight  line 
to  the  lowest  few  cross  section  points  and  extrapolating  to  zero  cross 
section. 

IV.  RESULTS  AND  DISCUSSION 

A total  of  ten  reactive  cross  section  curves,  a(v,  J;  E)  vs.  E, 
corresponding  to  the  initial  quantum  numbers  v = 0,  1 and  J = 0,  1,  2, 

3 and  4,  were  computed  for  the  reaction  O + H^Cv,  J) — *-OH  + H.  About 
300  trajectories  were  calculated  for  each  cross  section  point.  The 
coefficients  An(v,  J)  along  with  the  threshold  energies  E^(v,  J),  which 
when  substituted  into  Eq.  (5)  provide  a good  fit  to  the  calculated  cross 
sections,  are  listed  in  Table  II,  and  the  curves  are  displayed  graphically 
in  Figure  3. 

As  stated  above,  f(T)  is  the  probability  that  the  system  is  initially 

on  an  electronic  surface  that  will  allow  a reaction  to  occur.  The 

3 J + 

initial  state  of  the  reactants  is  0(  P)  + H,(  E ).  If  spin-orbit  forces 

are  neglected  this  state  is  ninefold  degenerate  and  correlates  with 

3 3 2 2 

both  E and  it  states  of  a linear  complex.  The  products  OH(  tt)  + H(  S) 

13  3 

correlate  with  tt  and  tt  states.  Thus,  the  tt  linear  complex,  which 
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Table  II 


Coefficients  and  Threshold  Energies  for  o(v,  J;  E)‘ 


-122. 460 
- 89.  974 

1040. 80 
680.39 

-2968.50 

-1740.80 

-157. 280 

1418. 30 

-4416.20 

- 88.781 

562. 99 

-1152.70 

- 71.766 

320. 91 

- 219.81 

- 52.  308 

193.32 

- 241. 11 

- 36.404 

106.59 

- 101.78 

- 47.  882 

137.45 

- 125.36 

- 39.  533 

77.44 

6.06 

- 47.327 

127. 71 

- 87.17 

8.03 
3.48 
3.  39 
3.  53 
3.  90 
5 


ENERGY  (kcal/mole) 


Figure  3.  A plot  of  the  reactive  cross  sections  ct(v,  J;  E).  The  numbers 
labeling  each  curve  are  the  rotational  quantum  number  J;  the 
set  of  curves  with  threshold  near  4.  0 correspond  to  the  vibra- 
tional quantum  number  v = 1 and  the  set  with  threshold  near 
8.  5 correspond  to  v = 0.  These  curves  are  a plot  of  the  analy- 
tic fit,  i.  e.  , Eq.  (5),  using  the  coefficients  listed  in  Table  II. 


has  six  states,  correlates  with  both  products  and  reactants.  This 
simple  analysis  indicates  that  six  of  the  nine  initial  states  correlate 
with  the  products  yielding  f(T)  = 6/9.  This  approximation  is  adequate 
at  high  temperature,  but  at  low  temperatures  the  splitting  of  the 

3 

initial  state  by  spin-orbit  forces  becomes  important.  The  P state 
splits  into  three  energy  levels.  The  lowest  level,  P2>  is  fivefold 
degenerate,  and  all  five  states  correlate  with  the  ground  state  of  the 
products.  The  next  level  is  3Pj,  which  is  158.5  cm"1  higher  in  energy 
and  is  threefold  degenerate.  One  of  these  states  correlates  directly 
with  the  ground  state  of  the  products,  and  we  assume  that  the  reaction 
can  occur  on  this  surface  but  cannot  occur  on  the  other  two  surfaces. 

3 i 

Finally,  the  PQ  state,  which  is  226.  5 cm"  above  the  lowest  state, 
does  not  lead  to  the  ground  state  of  the  products.  Assuming  this  sim- 
ple model  and  the  Boltzman  probability  of  occupying  the  three  levels, 
we  obtain  the  formula 


f(T)  = [5  + exp(-228/T)]/[5  + 3 exp(-228/T)  + exp(-326/T)]. 


This  function  varies  between  1 and  2/3  as  we  go  from  low  to  high 
temperatures.  Westenburg  and  de  Haas20  used  the  high  temperature 
value  of  2/3  in  their  absolute  rate  theory  calculations.  In  the  temper- 
ature range  of  interest  for  our  calculations,  298  - 1000 °K,  f(T)  varies 
from  0.  81  to  0.  71. 

The  calculated  rate  constants  k(v=0,  T)  and  k(v=l,  T)  are  dis- 
played graphically  in  an  Arrhenius  plot  in  Figure  4,  and  selected  points 
are  listed  in  Table  III.  [From  now  on  k(v=0,  T)  will  be  simply  desig- 
nated as  k,  and  k(v=l,  T)  will  be  designated  k^.]  In  the  temperature 
range  298-1000  K these  rates  are  fit  very  well  by  the  nonlinear 
Arrhenius  formulas 


k = 2.  81T  X 10"1*  exp(-4279/T) 
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Table  III 


t L.„ 


Computed  Results  at  Selected  Temperatures 


kt/k 


298 

4.  9 ( - 1 8)a 

2.7  (-  14)a 

5465 

4.4 

300 

5.4  (-18) 

2.  8 (-14) 

5173 

4.4 

320 

1.4  (-17) 

4.4  (-14) 

3107 

4.2 

350 

4.  8 (-17) 

7.8  (-14) 

1617 

4.0 

400 

2.  5 (-16) 

1.7  (-13) 

680 

3.  7 

450 

9.3  (-16) 

3.2  (-13) 

348 

3.  6 

500 

2.  7 (-15) 

5.5  (-13) 

204 

3.4 

700 

4. 3 (-14) 

2.3  (-12) 

52 

3.0 

1000 

3. 9 (-13) 

7.3  (-12) 

19 

2.8 

"I*  3 11 

aThe  units  of  k and  k are  cm  molecule  sec 


The  number  in 


parenthesis  is  the  power  of  10. 


-19- 





* 

I 


k1”  = 4.  65T  X 10"14  exp(-1868/T) 

3 -1  -1 

where  the  units  are  cm  molecule  sec 

| - 1 3 3 

Birely  et  al.'  have  set  an  upper  bound  on  k of  10  cm  molecule 

sec-1  at  300°K.  Our  value  at  this  temperature  is  2.  8 x 10"14  cm3 

molecule'*  sec  *.  This  upper  bound  is  3.6  times  larger  than  our 

computed  value.  Birely  has  also  set  an  upper  bound  on  the  ratio; 

k^/k  s 3.8  X 104.  Our  value  for  this  ratio  at  300  K is  k^/k  = 5.  2 X 103. 

13 

Birely's  ratio  was  calculated  using  the  extrapolated  Leeds  data  for  his 
value  of  k at  300°K,  whereas  our  results  are  computed  using  the  higher 
calculated  value  of  k at  this  temperature.  Values  of  our  computed  ratio 
at  selected  temperatures  are  listed  in  Table  III.  A good  fit  is  given  by 
the  formula 


kT/k  = 1.  65  exp(2411/T) 

The  branching  ratio  T was  calculated  for  the  case  in  which  is 
initially  in  the  v = 1 state.  This  quantity  is  by  definition  equal  to  ratio 
of  the  state -to -state  rate  constants,  k(v=l,  v'  = l;  T)/k(v=l,  v'=0;  T). 
These  rate  constants  were  evaluated  by  the  same  method  used  previously, 
i.  e.,  the  computed  state-to- state  cross  sections  aQ(v,  J,  v1;  E)  were  fit 
to  the  function  (5),  then  the  rates  were  evaluated  by  integrating  Eq.  (3) 
with  a(v,  J;  E)  replaced  by  oQ(v,  J,  v';  E).  Branching  ratio  values  at 
selected  temperatures  are  listed  in  Table  III.  These  points  can  be  fit 
quite  well  to  the  formula 

T = 2.  3 exp(l 96/T) 

The  value  of  the  branching  ratio  is  sensitive  to  the  method  used 
to  determine  the  final  vibration  quantum  number  v’ . The  method  used 
in  these  calculations  is  the  quasiclassical  trajectory  histogram  method 
in  which  v'  + 1/2  is  calculated'63  as  a continuous  variable  and  then 
rounded  to  the  nearest  integer.  This  could  probably  be  improved  on  by 
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using  the  recently  proposed  methods  of  Truhlar  and  Duff  and  of 

22 

Bowman  and  Leasure. 

9 t 

Experiments  are  in  progress  to  measure  the  rate  k and  the 
branching  ratio  T at  300  K.  The  results  are  still  being  analyzed,  so 
we  cannot  yet  give  a comparison  of  calculated  and  experimental  values. 

A great  number  of  experimental  measurements  of  k by  a variety 
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of  methods  and  over  a wide  temperature  range  have  been  carried  out.  ’ ' 

From  an  evaluation  of  selected  high  and  low  temperature  measurements, 
the  Leeds  group  (Baulch  et  aL),3have  recommended  the  expression 
k = 3 X 10" 14  Texp(-4480/T)  cm3  molecule"1  sec"1  in  the  temperature 
range  400-2000 °K.  Several  authors'2,23  have  recently  expressed  the 
opinion  that  this  formula  does  not  show  quite  enough  curvature;  in 
particular,  it  underestimates  the  rate  at  low  temperatures.  It  is 
plotted  as  the  dashed  curve  in  Figure  4.  Our  calculated  results  lie 
above  on  the  solid  curve  labeled  v = 0.  The  Campbell  and  Thrush  point 
is  also  plotted  in  this  figure.  At  500  K the  calculated  k is  larger  than 
the  Leeds  value  by  a factor  of  1.4  and  at  1000  K the  factor  is  1.  15. 

The  greatest  source  of  uncertainty  in  the  reliability  of  the  calcu- 
lated results  is  our  lack  of  detailed  information  about  the  potential 
surface.  Our  choice  of  the  simple  LEPS  function  is  a first  approxima- 
tion consistent  with  our  limited  knowledge  of  the  interatomic  potential. 

At  the  present  time  we  have  no  good  reason  for  choosing  any  other 
potential;  however,  it  is  useful  to  have  some  idea  of  the  effect  that 
certain  changes  in  the  shape  of  the  potential  can  have  on  the  computed 
rates. 

24 

The  extended  LEPS  function  is  a more  flexible  surface  with  three 

adjustable  (Sato)  parameters.  Two  of  these  parameters  will  be  equal 

due  to  the  symmetry  of  the  potential  with  respect  to  interchanging  the 

hydrogen  atoms.  This  leaves  two  quantities  to  adjust  which  allows  us, 

within  certain  limits,  to  control  both  the  position  and  height  of  the 

energy  barrier  independently.  The  effect  of  changing  the  barrier  posi- 
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tion  was  studied  by  Polanyi  et  al.  Their  results  show  that  for 
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Figure  4.  Rate  constants,  k(v,  T),  for  the  reaction  O + H2(v) — ►OH  + H 
where  v is  the  vibrational  quantum  number  of  H£.  The  dashed 
line  is  a plot  of  the  rate  k = 3.0  x 10"14  Texp(-4480/T),  recom- 
mended by  the  Leeds  group.*  Also  shown  are  the  experimental 
rate  (•)  of  Campbell  and  Thrush*1  and  the  experimental  upper 
bound  (X)  established  by  Birely  et  al.* 


endothermic  reactions,  moving  the  barrier  into  the  exit  valley  of  the 
potential  surface  increased  the  effectiveness  of  the  vibrational  energy 
in  promoting  reactions.  In  a preliminary  set  of  calculations  we  veri- 
fied these  general  conclusions  for  this  reaction.  The  Sato  parameters 
of  the  extended  L.EPS  surface  were  adjusted  (to  A_„  = 0.  14  and  A = 

0.  03025)  so  that  the  barrier  height  was  the  same  (12.  5 kcal/mole)  but 

° ° 

moved  further  into  the  exit  valley  (O-H  = 1,  07A  and  H-H  = 1.  03A  com- 
pared  to  O-H  = 1.12  and  H-H  = 0.  95  in  our  original  surface).  The 
results  showed  a small  increase  in  the  v = 0 reactive  cross  sections 
(about  10  percent  for  J = 1 at  10  kcal/mole)  and  a rather  substantial 
increase  in  the  v = 1 reactive  cross  section  (about  300  percent  for  J = 1 
at  6 kcal/mole)  along  with  a lowering  of  the  threshold  energy.  Several 
other  cross  sections  were  also  calculated  with  results  of  a similar 
order.  The  conclusion  is  that  both  k and  k^  increase  as  the  barrier 
moves  into  the  exit  valley  with  the  change  in  k^  large  compared  to  the 
change  in  k. 

Clearly,  further  work  is  required,  especially  with  regard  to  the 
potential  surface.  An  ab-initio  calculation  of  the  G + H^  reaction  sur- 
face  indicates  the  barrier  position  is  slightly  more  into  the  reactant 
channel  (O-H  = 1.  15A  and  H-H  = 0.  85A)  than  the  LEPS  surface  used 
here.  This  would  tend  to  diminish  the  vibrational  acceleration.  In 
order  to  more  accurately  determine  this  effect,  we  anticipate  carrying 
out  a new  set  of  trajectory  calculations  using  the  ab-initio  surface  when 
it  is  completed. 
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THE  IVAN  A.  GETTING  LABORATORIES 


The  Laboratory  Operation*  of  The  Aerospace  Corporation  i*  conducting 
experimental  and  theoretical  inve«tigation*  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scie-itific  devel- 
opments is  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  physics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemical  reactions 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronics  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metals. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radia- 
tion from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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